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ABSTRACT 


Context. Understanding the mechanisms that launch and shape powerful relativistic jets from supermassive black holes (SMBHs) in 
high-redshift active galactic nuclei (AGN) is crucial for probing the co-evolution of SMBHs and galaxies over cosmic time. 

Aims. We study the high-redshift (z = 3.396) blazar OH 471 to explore the jet launching mechanism in the early Universe. 

Methods. Using multi-frequency radio monitoring observations and high-resolution Very Long Baseline Interferometry imaging over 
three decades, we study the milliarcsecond structure and long-term variability of OH 471. 

Results. Spectral modelling of the radio flux densities reveals a synchrotron self-absorbed spectrum indicating strong magnetic fields 
within the compact core. By applying the flux freezing approximation, we estimate the magnetic flux carried by the jet and find that 
it reaches or exceeds theoretical predictions for jets powered by black hole spin energy via the Blandford-Znajek mechanism. This 
implies that OH 471 was in a magnetically arrested disk (MAD) state where the magnetic flux accumulated near the horizon regulates 
the accretion flow, allowing efficient extraction of black hole rotational energy. 

Conclusions. Our study demonstrates the dominance of MAD accretion in powering the prominent radio flares and relativistic jets 
observed in the radio-loud AGN OH 471 and statistical studies of large samples of high-redshift AGN will shed light on the role of 
MAD accretion in launching and accelerating the earliest relativistic jets. 
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1. Introduction 


Understanding the jet launching mechanism in high-redshift ac- 
tive galactic nuclei (AGN) is crucial for unravelling the complex 
interplay between supermassive black hole (SMBH) accretion 
and feedback processes in the early Universe (e.g., Di Matteo 
et al. 2005, 2008; Dubois et al. 2014). The magnetically arrested 
disk (MAD) model proposes that relativistic jets can be pow- 
ered efficiently when the magnetic flux threading the black hole 
reaches a critical threshold (Narayan et al. 2003). Previous stud- 
ies have analysed AGN jets at redshifts up to 2.37, with most 
below z ~ 2 (Zamaninasab et al. 2014), however, confirming or 
refuting the MAD paradigm in higher redshift blazars can pro- 
vide unique insights into the earliest epochs of SMBH accretion 
physics and jet formation mechanisms. 

In this paper, we study the radio source OH 471 (0642+449, 
J0646+4451), a high-redshift (z = 3.396) blazar that serves as 
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a unique laboratory for testing MAD models and probing jet 
physics in the early Universe. It was one of the first distant 
(z > 3) radio quasars identified in the 1970s (Gearhart et al. 
1972; Carswell & Strittmatter 1973; Crowther 1973). Early radio 
observations of OH 471 revealed a complex spectral structure 
with multiple synchrotron self-absorption (SSA) components on 
different physical scales peaking between 0.5—20 GHz (Gearhart 
et al. 1974; Jones & O’Dell 1974), indicating emission from 
several distinct regions in the relativistic jet. Very Long Base- 
line Interferometry (VLBI) observations established OH 471 as a 
blazar with a typical one-sided core-jet morphology on milliarc- 
second scales (Marscher & Shaffer 1980; Gurvits et al. 1992; 
Kellermann et al. 1998; Kovalev et al. 2005; Lister et al. 2013; 
Fey et al. 2015; Lobanov et al. 2015). 


Multi-frequency VLBI observations reveal a compact core 
region with low (2%) fractional polarization and transverse mag- 
netic fields (Lobanov et al. 2015). This contrasts with the sub- 
stantially higher fractional polarization in a jet component lo- 
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cated ~ 4 mas (~30 parsec) downstream the core, which dis- 
plays magnetic fields aligning with the jet. Furthermore, an ex- 
tremely high rotation measure (RM) of 9 x 10? rad m~? (redshift- 
corrected) is inferred in the inner jet just 0.5 mas from the 
core (Hovatta et al. 2012). This represents the largest redshift- 
corrected RM detected among similar high-redshift blazars. 

A key motivation of this study is to explore the jet launching 
mechanism in OH 471 by estimating the magnetic flux and test- 
ing MAD models. Evidence for sufficient magnetic flux would 
strongly support magnetically-driven relativistic jets in this high- 
redshift radio-loud quasar. 

A direct way to estimate the magnetic field strength is to fit 
the SSA spectrum of OH 471. The radio spectrum of OH 471 
displays dramatic evolution over decades, with the SSA turnover 
frequency shifting from ~1 GHz in early observations (Gearhart 
et al. 1974; Jones & O’ Dell 1974) to 10-20 GHz in recent mon- 
itoring data (Volvach et al. 2000; Sotnikova et al. 2021). This 
remarkable spectral variability likely implies substantial changes 
in the magnetic field strength, particle densities and opacity close 
to the SMBH. However, the spectrum has transitioned to an un- 
usually flat shape between 1-22 GHz since 2015, representing 
a substantial departure from the earlier SSA-dominated spec- 
trum. Such spectral flattening likely indicates the emergence of 
a newly generated jet component still embedded in the core. 

By connecting changes in the jet structure and kinematics on 
milli-arcsec (mas) scales to this long-term spectral evolution, we 
expect to reveal the key structural changes in the inner jet re- 
sponsible for the spectral change. Combining multi-epoch VLBI 
images over decades ! allow us to probe the jet kinematics. Iden- 
tifying newly ejected components associated with the spectral 
flattening since 2015 could reveal essential clues to the particle 
acceleration and magnetic field evolution near the central engine. 

In this paper, we provide an in-depth analysis of the jet prop- 
erties and magnetic flux using multi-frequency radio variability 
modelling and up to three decades of VLBI image data from 
1995 to 2022. The paper is organized as follows: Section 2 de- 
scribes the data; Section 3 presents the results; Section 4 pro- 
vides discussion; and Section 5 summarizes the conclusions. We 
assume a ACDM cosmology” with On = 0.27, Qa = 0.73, 
Ho = 70 kms"! Mpc"!. 


2. Observations and Data Analysis 


This study utilizes two main radio datasets to investigate the 
spectral evolution and jet structure of radio OH 471. 

The first dataset is comprised of radio flux density light 
curves monitored by the RATAN-600 radio telescope at the Spe- 
cial Astrophysical Observatory of Russian Academy of Sciences 
(Parijskij 1993). Continuous monitoring of OH 471 has been car- 
ried out since 1997 at six frequencies: 1.1, 2.3, 4.7, 7.7/8.2, 11.2, 
and 21.7/22.3 GHz (Mingaliev et al. 2014, 2017). The mean time 
interval is approximately 65 days. Five epochs of data points 
observed by the RT-32 Zelenchukskaya and Badary radio tele- 
scopes (Shuygina et al. 2019) at 5.05 and 8.63 GHz are also 
included. In total, the RATAN-600 and single-dish light curves 
(Figure 1) include over 105 epochs spanning 15 years from 2009 
to 2024 (Appendix A). This long-term multi-frequency monitor- 
ing provides excellent temporal sampling to quantify the evolu- 
tion of the radio spectrum. 


1 https:// astrogeo.org/ maintained by Leonid Petrov. 

> Using this cosmological model, 1 mas angular scale corresponds to 
a projected physical size of 7.635 pc at z = 3.396, and 1 mas yr! = 
109.5 c¢. 
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The second dataset consists of VLBI data at 2.3, 5, 8.4, and 
15 GHz observed between 1995 and 2022. The 2.3, 5 and 8.4 
GHz VLBI data are primarily obtained from the geodetic ob- 
servations archived in the Astrogeo database (e.g. Beasley et al. 
2002). The 15.4 GHz VLBI data are from the Monitoring Of Jets 
in Active galactic nuclei with VLBA Experiments (MOJAVE) 
programme ° (2). We used calibrated visibility data from these 
databases that had undergone initial calibration steps. We only 
performed a few iterations of self-calibration and made mode 
fitting to study the change of the jet structure over time. The de- 
tails of VLBI data processing are described in Appendix B. 

Combining multi-frequency light curves, VLBI imaging, and 
spectral modelling enables a robust investigation connecting the 
structural and spectral evolution of OH 471. 


3. Results 
3.1. Radio variability 


Figure | presents the flux density variations of OH 471 observed 
by VLBI (top panel) and RATAN-600 (bottom panel) over a 
wide range of radio frequencies. Two major flares are evident — 
around March 2003 and October 2008 — visible only at the higher 
frequencies of 8 and 15 GHz. The higher-frequency (> 4.7 GHz) 
flux densities gradually decrease after 2009, likely indicating a 
period of reduced activity. In contrast, the 2-GHz flux densities 
show a slow rise from 1997, peaks around December 2013, and 
declines thereafter. This behaviour likely reflects the shift in the 
SSA turnover frequency over time and the inherent lag in the 
low-frequency emission response to core flare activity. 

The RATAN-600 lightcurves start in October 2009, although 
the initial onset of the 2008 flare is not captured due to limited 
temporal coverage. Nevertheless, the decline from the peak flux 
density is well documented. Following this flare, the light curves 
exhibit significant complexity and irregularity, suggesting pos- 
sible overlapping events or a highly variable source structure. 
The 2008 flare declined until 2020, after which the flux densities 
showed a pronounced rise, particularly at the higher observed 
frequencies. This rise suggests the emergence of a new jet com- 
ponent. Continuous long-term monitoring is desirable to fully 
capture the evolution of this late-arising event, as it may analo- 
gously trace activities similar to that observed during the earlier 
prominent flares in the 2000s. Studying these prominent flares 
can help understand the possible connection between the flar- 
ing activity and the evolution of the jet structure in high-redshift 
blazars like OH 471. 


3.2. Jet structure and Kinematics 


Previous studies have revealed a jet extending 10 mas east of the 
core (Gurvits et al. 1992) and up to 30 mas at 1.4 GHz (MO- 
JAVE webpage) *, while higher frequency (v > 2.3 GHz) VLBI 
resolves more compact structure within ~2 mas (Marscher & 
Shaffer 1980; Kellermann et al. 1998; Kovalev et al. 2005; Lis- 
ter et al. 2013; Fey et al. 2015). Multi-epoch VLBI observations 
have revealed apparent superluminal speeds up to 8.5 c and con- 
strained key jet parameters including the Lorentz factor I = 5.4, 
viewing angle 6; = 0.8°, and intrinsic opening angle of ¢; = 0.3° 
(Pushkarev et al. 2009; Lister et al. 2013). These studies revealed 
both fast-moving and stationary jet features, indicating the coex- 


3 https: //www.cv.nrao.edu/MOJAVE/ 
4 https://www.cv.nrao. edu/MOJAVE/sourcepages/0642+449. 
shtml 
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istence of propagating shocks and recollimation sites or jet-ISM 
interactions. 

We analyzed over 50 epochs of 8.4-GHz and 2.3-GHz 
archival VLBI data spanning 1996-2022 from the Astrogeo 
database to qualify jet motions. The high resolution and long- 
time baseline provided by these observations enable robust con- 
straints on jet proper motions. We measure a proper motion of 
u = —0.0025 + 0.0085 mas yr! for component J1 located ~3.5 
mas from the core, and 0.040+0.003 mas yr“! (corresponding to 
4.4 + 0.3 c) for J2 at 0.5-1 mas from the core (Figure C.1). The 
stationary jet component J1 shows an extended morphology, and 
the polarized flux density along the southern edge of this fea- 
ture could suggest either an interaction between the jet and the 
ambient medium (Lobanov et al. 2015) or the presence of a he- 
lical magnetic field component roughly aligned with the jet. As 
the 8-GHz resolution differs from 15-GHz MOJAVE data, we do 
not directly compare our proper motion of J2 to literature values 
(Lister et al. 2013), which probe distinct jet structures. 

The close temporal association between the emergence of 
jet components and radio flares supports a direct connection be- 
tween the flaring activity and the ejection of a new jet component 
from the central engine. We found that J2 was ejected earlier than 
the monitored lightcurves (i.e. before 1996), based on backward 
extrapolation. The observed variability combines the ensemble 
core and jet emission. If each flare produces a component, but 
the VLBI resolution cannot distinguish them, the detected mo- 
tion represents overlapping components. This could explain why 
individual jet features do not always correspond to single flares. 
Alternatively, J2 may have decelerated during its outward mo- 
tion, causing an underestimation of its initial speed. Dividing the 
J2 data into pre- and post-2016 segments reveals faster motion 
(0.053 + 0.002 mas yr“) initially. Back-extrapolating this speed 
implies an ejection time around 2001, roughly consistent with 
the strongest observed flare during the past three decades. 


3.3. Radio spectra 


Long-term radio spectral monitoring of OH 471 has revealed sig- 
nificant variability (e.g., Volvach 2003). In the 1970s, the spec- 
trum showed a synchrotron self-absorption turnover around 1 
GHz and a flat spectrum at higher frequencies, possibly due 
to a combination of multiple SSA components (Gearhart et al. 
1974; Jones & O’Dell 1974). However, from the 1980s on- 
ward, this low-frequency SSA component faded away, while 
a new high-frequency SSA component emerged and became 
dominant. Since the 1990s, regular (1-22) GHz observations by 
RATAN-600 revealed strong turnover frequency evolution over 
time (Mingaliev et al. 2012; Sotnikova et al. 2021). The turnover 
was around 20 GHz in May 1998 (Volvach et al. 2000) and 
around 15 GHz in November 1998 (Dallacasa et al. 2000), drop- 
ping to 11.2 GHz by July 2002 (Tinti et al. 2005). Since 2012, it 
has remained variable around 17 GHz (Sotnikova et al. 2021). 

These long-term observations have revealed clear variabil- 
ity in both the flux density and peak frequency of the inverted 
spectrum of OH 471 over decades. To extract physical insights, 
we modelled each radio spectrum using a synchrotron self- 
absorption model (Appendix D), fitting for the turnover fre- 
quency Vm, flux density at the turnover Sm, the optically thick 
spectral index Qnick, and the optically thin spectral index ain 
(see a typical SSA spectrum in Figure 2). 

Figure 2 shows the remarkable evolution of the radio spec- 
trum from 2009 to 2024 based on multi-frequency monitoring 
data. The observations suggest at least two distinct SSA com- 
ponents: a low-frequency component (turnover < 5 GHz) cor- 


responding to the steady-state radio core, and a high-frequency 
flaring component (turnover > 5 GHz). This is consistent with 
conclusions from early studies (Gearhart et al. 1974). However, 
our data only allows fitting a single SSA model. The high- 
frequency component evolves faster and dominates the overall 
spectral changes. The new jet component associated with the 
2003 and/or 2008 flares gradually faded over 2009-2019, caus- 
ing Vm and ainick to decrease, until another new high-frequency 
component emerged after 2020, shifting the turnover back to 
higher frequencies. The joint decrease of vm and S m over 2009— 
2019 may imply a steady decline in the magnetic field strength 
and emitting particle densities in the core region. This is most 
straightforwardly explained by a reduction in the jet power from 
the central engine during this period. The combined changes in 
Qthick, Ym and Sm support the hypothesis of a major reconfigu- 
ration of the jet structure on sub-pc scales, with the formation 
of a new magnetized component. Currently, the flux densities of 
OH 471 are still in a rising phase. Closely monitoring this new 
component will help understand the details of jet production and 
evolution. 

Our approach of fitting a single SSA spectrum to the inte- 
grated radio flux densities has inherent limitations, as the ob- 
served flux contains contributions from both the core and jet 
components. Separating the core and jet reliably requires broader 
frequency coverage from VLBI imaging. By fitting a single spec- 
trum to the limited frequency data points, there are uncertainties 
introduced into the derived parameters like the magnetic field 
strength, though the core dominates the overall flux density in 
all epochs. We have attempted to fit the radio spectrum using 
two SSA components (i.e., a radio core and an inner jet), with 
the core being the dominant SSA component at high frequencies. 
Comparing the results with those from a single SSA component, 
we found that the flux density of the core decreases and the peak 
frequency increases. The combined effect of these changes leads 
to an increase in the revised Bssa, further supporting the paper’s 
conclusion that the jet magnetic flux in OH 471 is higher than 
in the MAD prediction. However, due to the current sparse num- 
ber of observational frequency points, we cannot accurately con- 
strain the SSA parameters. Although using the integrated flux 
density and a single SSA component to fit the radio spectrum is 
not perfect, it provides a way to constrain the source properties 
and magnetic fields in the vicinity of the black hole. These lim- 
itations should be kept in mind when interpreting the results of 
this study. Future multi-frequency VLBI studies can overcome 
this limitation by robustly separating the core and jet contribu- 
tions for more reliable spectral modelling. 


3.4. Brightness temperature of the VLBI core 


OH 471 stands out as one of the most luminous and beamed 
high-redshift blazars (e.g., O’ Sullivan et al. 2011). High bright- 
ness temperatures exceeding 10!* K observed for the OH 471 
core (Table B.1) are consistent with the expectations for a highly 
Doppler-boosted relativistic jet pointed roughly toward our line 
of sight (Section 3.2). We also find that the maximum brightness 
temperature Tp max Shows a systematic decreasing trend with in- 
creasing observation frequency: 


-— T, > 4x 10" K at 1.6 GHz from RadioAstron space VLBI 
observations (Lobanov et al. 2015); 

— Ty ~ 3.3 x 10!? K at 5 GHz from the VLBI Space Observa- 
tory Program observations (Dodson et al. 2008); 

— Ty = 1.2-4.3x 10! K at 15 GHz from ground-based VLBA 
observations (Kovalev et al. 2005); 
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- Ty = 1.1-1.6x10!! K at 86 GHz from the Global mm-VLBI 
Array (Lobanov et al. 2000; Lee et al. 2008). 


The decreasing trend in Tp,max of the OH 471 core with in- 
creasing observation frequency can be explained by a combi- 
nation of effects: synchrotron self-absorption causes the core to 
become more optically thin and extended at higher frequencies; 
the SSA turnover around 8 GHz (Section 3.1) suggests that emis- 
sions at frequencies below this peak are in the optically thick 
regime, contributing to higher observed temperatures; higher 
resolution available at higher frequencies may resolve substruc- 
ture and gradients in the core magnetic field, particle densi- 
ties, and relativistic beaming; different emission mechanisms can 
dominate across frequencies; unresolved sub-components likely 
contribute to the complex spectral behaviour. 


4. Discussion 


As one of the most distant blazars known, OH 471 provides a 
unique window into jet production and particle acceleration in 
the vicinity of an actively growing SMBH in the early Universe. 
Studying the variability, structural changes and magnetic field 
evolution of OH 471 is key to unravelling the physical processes 
driving powerful high-redshift jets. 

The magnetic flux paradigm proposes that jet power is funda- 
mentally determined by the accumulation of poloidal magnetic 
flux around the black hole (Sikora & Begelman 2013). Poloidal 
flux refers to the vertical magnetic field component that threads 
both the black hole and the inner accretion disk. In MADs, 
strong poloidal fields obstruct the inner accretion flow, chok- 
ing the accretion and creating a highly magnetized inner region 
(Tchekhovskoy et al. 2011). In MAD states, the strong vertical 
magnetic fields extract the spin energy of the black hole through 
the Blandford-Znajek process (Blandford & Znajek 1977), gen- 
erating jets with kinetic powers comparable to or exceeding the 
accretion power (Narayan et al. 2003). Without adequate mag- 
netic flux accumulation, AGN cannot achieve MAD conditions, 
resulting in weaker jets compared to accretion power. 

Observations support links between the magnetic flux and 
the jet power: e.g., MAD jets exhibit more ordered inner disk 
fields and higher core rotation measures (Zamaninasab et al. 
2014); intermittent jet activity implies fluctuating magnetic flux 
levels (Sikora & Begelman 2013); jet power correlates with radio 
core luminosity (Nemmen & Tchekhovskoy 2015). These man- 
ifest that creating powerful jets relies on accumulating strong 
poloidal magnetic flux, beyond just black hole spin. 

We can test if OH 471 has achieved the MAD state by com- 
paring its estimated jet magnetic flux to the predicted MAD flux, 
thus assessing the role of magnetic flux in powering the powerful 
high-redshift jet in OH 471. 

Models of jet launching propose that relativistic jets can 
be efficiently powered once the poloidal flux surrounding the 
black hole exceeds a critical threshold via the Blandford-Znajek 
mechanism. Simulations indicate that the MAD state occurs 
when magnetic flux reaches ~ 50 (Mr2c) (Tchekhovskoy 
et al. 2011), in which M is the accretion rate, rg the gravi- 
tational radius, and c the speed of light. The mass accretion 
rate is M = Lice /(ne?), where Lacc is the accretion power, 7 
the efficiency (typically 0.1-0.4). The black hole mass (Mgy) 
sets the gravitational radius via rg = GMpy/c?, where G is 
the gravitational constant. Combining these quantities, the pre- 


dicted MAD magnetic flux is yap ~ 50 4| Laer (He =24x 
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nosity Lace is calculated using the bolometric luminosity Leo; and 
an accretion efficiency, giving Lace = Peek The bolometric lu- 
minosity was estimated from the optical C IV line luminosity. 
OH 471 presents strong optical emissions, which can be used to 
estimate Lace ~ 10LELR = 88.2Lc 1 ¥ 6.5 x 10% erg s7! (Liu 
et al. 2006; Francis et al. 1991; Ghisellini & Tavecchio 2008). 
With 7 ~ 0.4 (Zamaninasab et al. 2014), the estimated MAD 


flux is map © 1.7 x 10°, 484] [G cm’]. 


While we lack direct observations of the jet magnetic flux 
in OH 471, we can estimate it from synchrotron self-absorption 
magnetic field strength (Bssq) derived from the spectrum fitting 
and core size measurement. We estimate the SSA magnetic field 
strength following Marscher (1983) with details presented in Ap- 
pendix E. Bssa depends on turnover frequency, peak flux den- 
sity, core size, and the Lorentz factor (Appendix D). The core 
size is derived from VLBI model fitting (Appendix B). Follow- 
ing Chamani et al. (2021), Oje = 0.8X 107 f(ax)(1 +0) PTE] 
[G cm?], where o is the jet magnetisation parameter, f(a*) is 
a dimensionless parameter related to the black hole spin, B is 
the magnetic field strength at 1 parsec from the black hole. 


f(a*) = [1+ V1 — a?]/ais a function of the black hole spin a. For 
a rapidly rotating black hole, which is expected for powerful jet 
sources such as blazars, f(a*) ~ 1 (Appendix F.2). The magnetic 
field strength at a distance of 1 pc from the black hole, Bipc, was 
estimated using the relation Bipe = Beore X (Acore/ Ipc), where 
heore 18 the distance between the VLBI core and the central en- 
gine, Beore is the core magnetic field strength derived from the 
SSA fitting, and £ is the power-law index describing the mag- 
netic field decay with distance, as expected for a conical jet ge- 
ometry (?). In this study, we adopt 8 = 1, which corresponds 
to the magnetic field being inversely proportional to the distance 
from the black hole. Given the conical jet geometry, the angu- 
lar size of a spherical homogeneous synchrotron source at the 
turnover frequency (ym) is approximately equal to the angular 
size observed at another frequency 6(v) multiplied by the ratio of 
the frequencies. Further details on the calculation of Bıpc can be 
found in Appendix E. 


To independently estimate the magnetic field strength in the 
VLBI core region, we performed a joint analysis of the simulta- 
neous 8.4 and 15.4 GHz images on 2011 March 5 to determine 
the core shift. We found a typical core shift of (0.09 + 0.01) mas 
between these frequencies. Using the method given by ?, assum- 
ing equipartition between the particle and magnetic field energy 
densities, and assuming a viewing angle of 0.25°, we estimate a 
magnetic field strength at a distance of 1 pc from the black hole 
of Bipe = (4.2 + 0.8) G. This is slightly lower than the value 
of (10 + 2) G derived from the SSA fitting, but the two values 
agree to within 30. The core shift estimate thus provides an in- 
dependent confirmation of the high magnetic field strength in 
the OH 471 jet, although both methods rely on several simplify- 
ing assumptions, so the true uncertainties are likely to be larger 
than the formal errors reported here. In addition, we have only 
one epoch of quasi-simultaneous data, so the core-shift calcu- 
lation here represents a rough assessment of the magnetic field 
strength derived from the SSA spectrum fit. Higher-resolution 
multi-frequency VLBI observations are desirable to refine these 
magnetic field estimates and to probe the magnetic field structure 
of the jet in more detail. 

The calculated Bssa, Bipc, and magnetic flux Pjes are shown 
in Figure 4. Bssa does not exhibit systematic variations over 
time. Bssa depends on multiple parameters (Bssa œ ve S zô) 
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that vary differently over time. These combined effects essen- 
tially cancel out, leading to a relatively stable Bssq. Physically, 
Bssa represents the magnetic field strength right at the SSA pho- 
tosphere. This photosphere likely moves inward/outward as the 
jet conditions evolve, sampling different field strengths across a 
range of radii. In contrast, Bıpc approximates the field strength at 
a fixed 1 parsec distance. This single location maintains a higher 
field when the jet is actively flaring with fresh particle accelera- 
tion. As the flare fades over time, the field decreases at that fixed 
point. As Bj). is more heavily dependent on vm through the ve 
scaling, Bıpc follows a similar evolution trend with vm. 

The inferred Ọja exceeds the predicted Oman in 2009-2012, 
when the source was in a declining phase after the large 2008 
flare and the newly generated jet had not yet dimmed. From 
2015-2020, Pj, ~ map, indicating magnetic flux accumula- 
tion between flares. The period of 2021 December — 2024 Febru- 
ary represents the rising phase of a new flare whose peak is yet 
to be determined, and in this period Ọje is significantly higher 
than ®yap, indicating that the magnetic flux has now exceeded 
the critical threshold needed to trigger a new jet ejection episode. 
The strong magnetic flux explains how OH 471 is able to launch 
such a powerful jet. 


The temporal evolution of Ọje relative to Dyap implies a 
cyclical flux accumulation process, where sufficient magnetic 
flux builds up to launch a new jet component, followed by a de- 
clining phase as the new jet expands and fades. Monitoring the 
magnetic flux over multiple flaring cycles can further test this 
magnetic paradigm for powering recurring jets. 

To maintain the MAD state for long periods of time, signifi- 
cant magnetic flux is likely to have accumulated early on near the 
event horizon of the black hole. This type of accretion is thought 
to fuel the rapid growth of SMBHs at high redshifts. This sce- 
nario suggests a link between the jet’s current activity and the 
AGN’s early accretion history. The large reservoir of magnetic 
flux from past accretion builds up over timescales much longer 
than the jet variability, enabling the AGN to periodically reach 
MAD conditions and drive powerful jets intermittently. Despite 
fluctuations in magnetic flux, the stored flux near the black hole 
ensures continued jet activity following accretion episodes. 

In Figure 4, we plot the measured jet magnetic flux Ọjet ver- 


sus Le Mpu for OH 471 compared to the sample from Zamani- 
nasab et al. (2014). Their sample consists of 76 radio-loud AGN, 
including 68 blazars and 8 nearby radio galaxies. Zamaninasab 
et al. (2014) used published measurements of the core-shift ef- 
fect to calculate the jet magnetic flux for each source. Since jet 
varies with time and Mgy is not determined for OH 471, we show 
the range that spans the Ọje variations and Mgy = 108710 Mo. 
The (Ọja) value derived for OH 471 broadly aligns with the 
previously established relationship between jet magnetic flux 
and the combined quantity (L! Mpp) for powerful jets from 
z = (0.07 — 2.37). While our calculation involves uncertain- 
ties discussed before in this section, the overall consistency with 
lower-redshift AGN jets suggests the magnetic flux paradigm re- 
mains applicable to powerful jets at higher redshifts. The agree- 
ment between OH 471 and the relation for lower-z sources hints 
at the potential universality of magnetically-driven launching 
processes of powerful jets from the local Universe to earlier 
cosmic epochs. Further studies of additional high-redshift jet- 
ted AGN can solidify the connection between magnetic flux and 
jet power from the early AGN population to their modern coun- 
terparts. 


The magnetic field estimates presented in this paper are 
based on the SSA spectrum fitting of the entire source. How- 


ever, this approach has some limitations due to the mixing of the 
core and inner jet components at lower frequencies (e.g. 2.3 and 
5 GHz), which can introduce uncertainties in the derived mag- 
netic field strengths. To obtain more accurate magnetic field esti- 
mates and to probe the magnetic field structure of the jet in detail, 
high-resolution multi-frequency VLBI observations are needed. 
These observations will allow us to more precisely determine the 
core shift between different frequencies, providing an indepen- 
dent constraint on the magnetic field strength. By comparing the 
magnetic field estimates derived from the core shift analysis with 
those obtained from the SSA spectrum fitting, we will be able to 
assess the robustness of our results and better constrain the un- 
certainties. Furthermore, the high-resolution images will allow 
us to study the evolution of the magnetic field along the jet, pro- 
viding insights into the collimation and acceleration processes of 
the jet. 


5. Conclusion and Summary 


As a radio-loud quasar at z > 3, OH 471 provides a unique 
window into jet formation and evolution in the early Universe. 
It shows multiple observational signatures that are consistent 
with the predictions of a magnetically arrested disk based on 
the magnetic flux paradigm for powering relativistic jets — a sys- 
tem in which strong magnetic flux accumulation enables efficient 
extraction of black hole spin energy to power relativistic jets. 
These indications include: extremely high brightness tempera- 
tures exceeding 10! K measured with VLBI on parsec scales 
(Lobanov et al. 2000), implying efficient jet formation from the 
core region; large rotation measures detected in the inner jet 
(Hovatta et al. 2012; O’Sullivan et al. 2011), which could arise 
from ordered magnetic field components and/or dense magne- 
tised plasma screens local to the jet formation region; long-term 
variability in the radio flux density and spectral shape (Lister 
et al. 2019; Sotnikova et al. 2021, and the present paper), consis- 
tent with the recurrent ejection of new synchrotron self-absorbed 
components seen in many AGN jets including OH 471, which 
potentially arise from episodic accumulation and release of mag- 
netic flux driving transient MAD episodes near the black hole; 
and an estimated jet magnetic flux that surpasses theoretical pre- 
dictions for the critical MAD threshold needed to power a rel- 
ativistic jet (the present paper). All these observations support 
a scenario where cyclic magnetic flux accumulation likely en- 
ables this high-redshift blazar to periodically achieve magneti- 
cally arrested states that give rise to energetic jet ejections from 
the vicinity of the central SMBH. 

This finding provides important insights into existing mod- 
els of the earliest growth of supermassive black holes. Power- 
ful flares and extremely relativistic jets observed in high-redshift 
AGN likely require achieving MAD conditions early, allowing 
time to accumulate magnetic flux. We expect other high-z jetted 
AGN to follow this magnetic paradigm, with the most powerful 
systems already reaching MAD levels. Studying a large statis- 
tical sample for correlations between jet power, accretion rate, 
and magnetic flux can firmly establish the importance of mag- 
netic flux accumulation in generating the earliest AGN jets. Fu- 
ture multi-band observations are crucial to understanding if the 
jet dynamics in OH 471 differ from similar lower-z objects, po- 
tentially revealing changes in jet production and collimation over 
cosmic time. 

While the MAD model provides rich qualitative predictions 
for the accretion and jet dynamics around black holes, it does 
not provide universal quantitative relations directly linking mag- 
netic flux accumulation and jet cycles to system parameters. The 
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complex interplay between black hole mass, accretion rate, disk 
properties, and GRMHD processes makes accurate quantitative 
modelling a challenge. Nevertheless, the MAD framework pro- 
vides important insights into the physical mechanisms likely to 
drive the episodic powerful jets observed from sources such as 
OH 471. Despite the limitations of quantitative modelling of in- 
dividual objects, multi-wavelength observations may help to dis- 
tinguish MAD signatures. 
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Fig. 1. Radio light curves of OH 471 made with the VLBI (top panel) , RATAN-600 and RT-32 (bottom panel) observations. 


{— prior 


— mcmc-fitting 


Flux density [Jy] 


DATE 2017.10. 1 Radio spectrum evolution 


2009-10-15 
2009-10-16 
2009-11-15 
2010-08-19 
2011-01-15 
2011-10-13 
2011-12-25 
—— 2017-04-10 
—— 2017-10-01 


Flux density [Jy] 


Vm = 6.47 + 0.62 GHz 


Sm = 2.26 + 0.08 Jy 


2017-12-18 
2018-04-29 
2018-07-13 
2019-10-09 
2021-12-07 
2022-11-05 
2023-11-07 
2024-02-20 


Frequency [GHz] 


Frequency [GHz] 


Fig. 2. Examples of Synchrotron Self-Absorption model fit (left panel) and the evolution of the fitted spectra (right panel). In the left panel, red, 
gray, and black lines denote the best mcmc-fitting result, 5000 sets of models randomly selected in the parameter space, and the prior parameters, 
respectively. The right panel collected the fitted results of 17 epochs from 2009 to 2024, with darker lines representing more recent observations. 
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Fig. 3. Variations of the jet physical conditions with time. The upper four panels illustrate the time-dependent evolution of the synchrotron self- 
absorption spectrum parameters. The lower three panels display the calculated magnetic field strength (Bssq and Bıpc) and magnetic flux, providing 
insights into the magnetic environment change over time. The horizontal line in the bottom panel denotes the predicted MAD magnetic flux. 
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Fig. 4. Measured magnetic flux versus Li? Mpu. The black circles are taken from Zamaninasab et al. (2014), consisting of 76 radio-loud AGN. 
The parameter space for the OH 471 data is shown as a shaded area, and the data points on three representative epochs are shown as large circles. 
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Table A.1. RATAN-600 continuum radiometer parameters: the central 
frequency fo, the bandwidth Afo, the detection limit for point sources 
per transit AF. FWHMraxpec is the angular resolution along RA and 
Dec calculated for the average angles. 


fo Afo AF FWHMpaxpee. 
(GHz) (GHz) (mJy beam!) 
22.3 2.5 50 0/17 x 1/6 
11.2 1.4 15 0:34 x 3/2 
8.2 1.0 10 0:47 x 44 
4.7 0.6 8 0/81 x76 
2:2) 0.08 40 1:771 x 15/8 
1.25 0.08 200 3/07 x 27/2 


Appendix A: RATAN-600 and RT-32 light curves 


The RATAN-600 radio telescope, a 600-meter ring-shaped 
multi-element antenna operating in transit mode (Korolkov & 
Pariiskii 1979; Parijskij 1993). This design allows for instanta- 
neous broadband spectral observations from 1—22 GHz within 3— 
5 minutes when a source transits across the focal line where the 
receivers are located. The angular resolution, measured by the 
full width at half maximum (FWHM), depends on the antenna 
elevation angle, and the resolution along declination FWHM pec 
is three to five times worse than that along the right ascension 
resolution FWHMka. The angular resolutions along RA and Dec 
calculated for the average angles are presented in Table A.1 for 
the six observing frequencies. Under favourable conditions and 
average elevation angles (Dec ~ 0°), the detection limit reaches 
around 5 mJy at 4.7 GHz for a RATAN-600 single sector for an 
integration time of ~ 3s. 

The measurements were processed using an automated data 
reduction system (Udovitskiy et al. 2016; Tsybulev et al. 2018) 
and the Flexible Astronomical Data Processing System (FADPS) 
standard package modules (Verkhodanov 1997) for the broad- 
band RATAN-600 continuum radiometers. Flux density calibra- 
tion relied on seven well-established secondary standards which 
are selected due to their high radio brightness and stable flux 
densities: 3C 48 (An et al. 2010), 3C 147, 3C 286 (An et al. 
2017), NGC 7027, DR 21, 3C 295 and 3C 309.1. The flux den- 
sity scales were calculated based on Baars et al. (1977) and Per- 
ley & Butler (2013, 2017). The measurements of the calibrators 
were corrected for the angular size and linear polarization based 
on the data from Ott et al. (1994) and Tabara & Inoue (1980). 

The total flux density uncertainty accounts for both the 
RATAN-600 calibration curve error, and the random antenna 
temperature measurement error (Udovitskiy et al. 2016), calcu- 
lated as: 


(=) = (zi) GS) 
peas =j + | —— 3 
S v g(h) Tanty 
where ors is the total flux density standard error; S,, the flux 
density at a frequency y; oc, the standard calibration curve error, 
which is about 1-2% and 2-5% at 4.7 and 8.2 GHz respectively; 
g(h), the elevation angle calibration function; om, the standard 
error of the antenna temperature measurement; and Tant is an 
antenna temperature. The systematic uncertainty of the absolute 
flux density scale (83-10% at 1-22 GHz) is not included in the 
total flux error. 

The RT-32 radio telescopes of the Institute of Applied As- 
tronomy of the Russian Academy of Sciences (IAA RAS) at 
the Zelenchukskaya (Zc) and Badary (Bd) stations were used to 
measure the flux densities at frequencies of 5.05 and 8.63 GHz, 


(A.1) 
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with the bandwidth A fo>=900 MHz for both receivers. The obser- 
vations were conducted in the drift scan mode. One scan lasted 
~] minute at 8.63 GHz and ~1.5 minutes at 5.05 GHz with 1s 
registration time, the FWHM=3’9 and 7/0 at 8.63 and 5.05 GHz 
respectively; flux density limit reaching of about 20 mJy per scan 
for both frequencies under optimal observation conditions. To in- 
crease sensitivity, the scanning cycle was repeated multiple times 
to form a continuous observing set. 

The observed data were processed using the original pro- 
gram package CV (Kharinov & Yablokova 2012) and the 
Database of Radiometric Observations. Scans corrupted by 
weather or interference were filtered out. The remaining scans 
were averaged and fitted with a Gaussian curve after baseline 
subtraction using a parabolic approximation. The antenna tem- 
perature and its error were estimated from the Gaussian analy- 
sis of the averaged scan. The reference signal error of the noise 
generator is less than 1% and is also included in the result. The 
flux density scales were calculated similarly to the RATAN-600 
observations (Baars et al. 1977; Perley & Butler 2013). 3C 48, 
3C 147, 3C 295 and 3C 309.1 were used as the reference sources. 

Figure 1 bottom panel shows the RATAN-600 lightcurves of 
OH 471. Several data points observed by RT-32 at 8 and 5 GHz 
after September 2022 are added in the lightcurve. Data points 
that significantly diverge from the expected evolutionary trend, 
likely due to observational errors, are excluded from the analysis. 


Appendix B: VLBI Observations and results 


The visibility data have been calibrated using standard proce- 
dures in the PIMA (Petrov et al. 2011) and AIPS software pack- 
age (Greisen 2003). We only performed hybrid imaging inter- 
actively using self-calibration and deconvolution cycles in the 
DIFMAP software package (Shepherd 1997). The final images 
were made using natural weighting after several rounds of phase- 
only self-calibration followed by 1-2 amplitude self-calibration 
runs. Data severely affected by Radio Frequency Interference 
(RFI), insufficient integration time or sparse (u, v) coverage, 
which could not yield high-quality images, were excluded from 
our analysis. The VLBI data analysis, as well as other data anal- 
yses in this research, were used vlbi-pipeline ° deployed on the 
computing platform of the China SKA Regional Centre (An et al. 
2019, 2022). 

Following the imaging process, we fit circular Gaussian 
models to the self-calibrated visibility data for each epoch us- 
ing the MODELFIT task in DIFMAP. This enabled a quantita- 
tive characterization of the milliarcsecond (mas)-scale core and 
jet components, including their flux densities, sizes, and relative 
positions over time. Elliptical or circular Gaussians were fitted to 
the self-calibrated visibility data. Typically 1-2 jet components 
in addition to the core are adequately used to model the visibility 
data at each epoch. The fitted 8.4-GHz Gaussian parameters are 
summarized in Table B.1. The uncertainties of the fitted parame- 
ters were estimated following (Li et al. 2020). By registering the 
jet component positions relative to the core, we measured proper 
motions to analyze jet kinematics. 

VLBI observations spanning decades have provided detailed 
insights into the parsec-scale structure of the blazar OH 471. 
Early 18 cm (1.67 GHz) observations revealed a compact 2.2 
mas core along with hints of more extended emission that could 
account for the flux density offset between VLBI and single-dish 
measurements (Marscher & Shaffer 1980). Subsequent VLBI 
monitoring firmly established OH 471’s characteristic core-jet 
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morphology on milliarcsecond (mas) scales (Gurvits et al. 1992; 
Kellermann et al. 1998; Kovalev et al. 2005), with extremely 
high brightness temperatures exceeding 10! K inferred from 
1.6 GHz space VLBI observations (Lobanov et al. 2015). High- 
frequency (15-43 GHz) ground-based and 1.6 GHz space VLBI 
images resolve the core into two components separated by 0.76 
mas along a position angle of ~ 81° east (Lister et al. 2013; Fey 
et al. 2015; Lobanov et al. 2015). While lower-frequency global 
and space VLBI observations detected an extended jet compo- 
nent stretching up to 30 mas (Gurvits et al. 1992; Lobanov et al. 
2000), higher-frequency VLBA images have failed to detect this 
distant emission, potentially due to resolution effects or a steep- 
ening spectral index. 

Complementing these findings, the VLBI images presented 
in Figure B.1 provide finer scale views of the core-jet morphol- 
ogy within 2 mas across a range of frequencies, showing similar 
core-jet structures at different frequencies. Figure | shows the ra- 
dio lightcurves derived from VLBI observations, which account 
for the total flux densities integrated over the entire core-jet re- 
gion. 


Appendix C: Jet kinematics 


Figure C.1 displays the proper motions of components J1 and 
J2. J1 does not show significant proper motion. In contrast, the 
inner jet J2 exhibits an outward motion of 0.040+0.003 mas yr”! 
(4.4 + 0.3 c). As J1 was not detected proper motion based on 2.3 
GHz VLBI data, their information is not presented in Table B.1. 

To evaluate the relativistic beaming effect inside the target 
OH 471, we estimated the jet kinematics properties in different 
epochs. We determined the core brightness temperature using the 
equation (e.g., Condon et al. 1982): 


Tpobs = 1.22 x 10'7(1 +2) 


¥ Prnailenin , 

where the flux density S is in the unit of Jy, the observing fre- 
quency v in GHz, and @ is the elliptical/circular Gaussian com- 
ponent angular diameter (full-width half-maximum, FWHM) in 
mas. The brightness temperature results at 8 GHz are listed in 
Table B.1. From the measurement of brightness temperature, we 


could estimate the Doppler factor using 6 = Pe , where Tp int 
represents the intrinsic brightness temperature. Assuming the en- 
ergy equipartition between radiation particles and the magnetic 
field, Tint ~ 5 Xx 10!°K (Readhead 1994). The corresponding 
Doppler factors are 25.0 < 6 < 49.9. 

The features of the jet, bulk Lorentz factor, I, and viewing 
angle to the line of sight, 0, can be calculated from the following 
formulae (Urry & Padovani 1995): 


B+eP+1 
[= -—*+__., 
26 
2Ba 
pip. 
Bi+e-1 


The £4 in the formulae is jet proper motion speed. We measured 
from VLBI long-term monitoring at X-band and constrain this 
speed to be 8, = 4.4 +0.3c. Adopting this result and the Doppler 
factor in different epochs, the bulk Lorentz factor and viewing 
angle are 12.4 ¢ T g 25.2 and 0.2° < 0 < 0.9°. These param- 
eters are used in calculating the magnetic field strength for the 
corresponding RATAN observation epochs. 


Appendix D: Synchrotron self-absorption spectrum 
model fitting 


The synchrotron self-absorption (SSA) spectrum is character- 
ized by a distinctive shape: it rises steeply at low frequencies 
where the emission is self-absorbed, reaches a peak where the 
optical depth is approximately unity, and then falls off at higher 
frequencies where the emission becomes optically thin. The con- 
ventional SSA spectrum can be described by the following equa- 
tion: 


mene ete) 


where S, represents the flux density at frequency v; Sm is the 
peak flux density at the turnover frequency; vm is the turnover 
frequency; a is the spectral index in the optically thin part of the 


(D.1) 


spectrum. The term f=)" describes the optically thick part of 
the spectrum where v < vm, and the emission rises steeply with 
increasing frequency. The term in square brackets accounts for 
the transition from optically thick to thin emission, and the ex- 
ponential term represents the attenuation of the emission due to 
self-absorption. The modelled vm constrains the electron density 
and magnetic field strength at the emission site. The peak flux 
S m indicates total radiative power. The index ainin provides clues 
about particle acceleration and cooling mechanisms. The exact 
shape of the SSA spectrum can vary depending on the details of 
the electron energy distribution and the geometry of the emitting 
region. The spectral index œ can also provide information about 
the energy distribution of the radiating electrons, which is often 
assumed to follow a power law. When fitting an SSA spectrum 
to observational data, the values of Sm, Ym, and œ are typically 
free parameters that are adjusted to obtain the best fit to the data. 
This fit can then be used to infer the physical properties of the 
source, such as the magnetic field strength and electron density, 
as previously discussed. 

Tiirler et al. (2000) derived an SSA spectrum equation (their 
equation 41): 


y Qthick 1 = exp (-tm (zy) 
1 — exp(-Tm) 


5,=Sn( (D.2) 


Ym 


Bahin 
Jorm 1} 

The main difference between this equation D.2 and the con- 
ventional SSA spectrum equation D.1 is the introduction of Tm as 
an approximation of the optical depth at the turnover frequency 
Vm, and as a function of the spectral indices @hick and Qthin. In 
the conventional SSA spectrum (Eq. D.1), the transition between 
the optically thick and thin regimes is often modelled with a 
simpler power-law dependence on frequency without explicitly 
accounting for the optical depth in this way. The equation D.2 
represents a more sophisticated approach to modelling the SSA 
spectrum by accounting for the frequency dependence of optical 
depth, which can lead to a more accurate characterization of the 
physical conditions within the source, especially when the spec- 
trum does not follow a simple power-law behaviour across the 
turnover point. This change is made more reasonable for AGN 
by taking into account an inhomogeneous jet consisting of mul- 
tiple homogeneous sub-regions, each with its own turnover fre- 
quency, resulting in a flatter optically thick spectral index than 
what would be expected from a homogeneous source. This in- 
homogeneity allows for a smoother transition between the op- 


where, Tm © 3 ( 
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Fig. B.1. Representative VLBI images of OH 471 across frequencies: This figure displays VLBI images at close epochs from right to left at 
frequencies of 15.4, 8.7, 2.3, and 1.4 GHz. Image details: 15.4 GHz with beam FWHM of 0.9 mas x 0.7 mas, PA = 6.9°, a peak flux density of 
3.04 Jy beam™! and rms noise of 0.98 mJy beam™!; 8.7 GHz with beam FWHM of 1.0 mas x 0.8 mas, PA = 0.4°, a peak flux density of 2.57 Jy 
beam! and rms noise of 0.72 mJy beam™!; 2.3 GHz with beam FWHM of 4.1 mas x 2.9 mas, PA=2.4°, a peak flux density of 1.474 Jy beam™! 
and rms noise of 0.85 mJy beam7!; and 1.4 GHz with beam FWHM of 10.9 mas x 7.6 mas, PA=25.6°, a peak flux density of 0.83 Jy beam™! and 
rms noise of 0.15 mJy beam~!. The lowest contour is set at three times the rms noise level, with contours increasing in a step of 2. 
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Fig. B.2. Changes of core size of OH 471 over time. A dashed line represents the average core size (0.30 + 0.11 mas) from 8-GHz VLBI 
observations. Notably, three data points showing a core size greater than 0.4 mas are linked to observations with limited (u, v) coverage, where 


only the core, lacking any jet structure, was detected. This observational limitation influences the measured core sizes. 


tically thick and thin regimes of the spectrum. The assump- 
tions of the model include a random orientation of the magnetic 
field within the jet, an electron energy distribution that follows a 
power law, and the application of general synchrotron radiation 
theory that allows the summing of emissions from various sub- 
regions. Moreover, relativistic effects such as Doppler boosting 
are assumed to be already included in the parameters, reflecting 
the complex motions within relativistic jets such as those from 
AGN. In addition, the use of the Tm approximation simplifies the 
determination of the optical depth at the turnover frequency. The 
inclusion of the turnover flux density S m as a separate parameter 
allows for a direct fit to the observed peak of the spectrum, so 
that the model fitting accurately captures the peak of the SSA 
spectrum. 


In this work, we derive the SSA model parameters by apply- 
ing the Bayesian Markov chain Monte Carlo (MCMC) algorithm 
with the Python package emcee (Foreman-Mackey et al. 2013). 
The input prior parameter set is fitted from Linear Regression 
with the SSA model equation. Then we randomly select 5000 
sets of parameters in the parameter space. The boundaries of the 
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four parameter is: 0 < Sm < 5,0 < vm < 20,0 < thick < 3, 
and —3 < @thin < 0. We then evaluate the fitting results with cor- 
ner figures and mean square error (MSE). We select the spectrum 
with a critical peaked profile for the magnetic field measurement. 


Appendix E: SSA magnetic field strength 


Equation D.2 not only provides a phenomenological fit but can 
also be used to derive a direct physical interpretation. The mag- 
netic field strength (B) can be estimated using the formula de- 
rived from the theory of synchrotron emission, which incorpo- 
rates the turnover frequency (vm), the source size, and the as- 
sumption of equipartition of energy between the magnetic field 
and the particles (Marscher 1983): 


4 ~2 
Bssa = 107b(a (2) (E) (3) (; —) 


where z is the redshift of the source; Bssa represents the mag- 
netic field strength in Gauss (G); b(œ) is a dimensionless param- 
eter depending on the optically thin spectral index and ranges 


(E.1) 
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Table B.1. Parameters of the fitted Gaussian model components 


Epoch Freq Comp S R OFWHM Ty 
(GHz) (mJy) (mas) (mas xmas) (10!!K) 

(1) (2) (3) (4) (5) (6) (7) 

2008/04/02 8.65 C 2.74 + 0.14 0.36 x 0.10 55.8 
J2 0.46 + 0.02 0.36 0.46 x 0.46 

2008/07/09 8.65 C 2.93 + 0.15 0.37 x 0.18 31.9 
J2 0.34 +0.02 0.42 0.51 x 0.51 

2009/05/13 8.65 C 2.76 + 0.14 0.37 x 0.18 29.3 
J2 0.28 +0.01 0.47 0.55 x 0.55 

2012/06/27 8.64 C 2.22 + 0.11 0.35 x 0.35 13.1 
J2 0.52 +0.03 0.62 0.43 x 0.43 

2013/12/11 8.64 C 1.86 + 0.09 0.54 x 0.54 4.5 
J2 0.40 +0.02 0.62 0.67 x 0.67 

2014/12/20 8.67 C 1.64 + 0.08 0.54 x 0.15 14.9 
J2 0.46 + 0.02 0.84 0.43 x 0.43 

2015/12/08 7.62 C 1.45 + 0.07 0.33 x 0.33 9.4 
J2 0.67+0.03 0.82 0.41 x 0.41 

2016/06/25 7.62 C 1.37 + 0.07 0.29 x 0.29 11.9 
J2 0.80 +0.04 0.78 0.51 x 0.51 

2017/01/16 8.65 C 1.42 + 0.07 0.51 x 0.17 11.8 
J2 0.60 +0.03 0.84 0.55 x 0.55 

2017/01/21 8.67 C 1.43 + 0.07 0.53 x 0.13 15.2 
J2 0.42 +0.02 0.90 0.35 x 0.35 

2017/05/27 8.67 C 1.41 + 0.07 0.45 x 0.45 5.0 
J2 0.69 + 0.03 0.83 0.57 x 0.57 

2018/05/07 8.65 C 1.34 + 0.07 0.64 x 0.16 9.7 
J2 0.61 +0.03 0.87 0.58 x 0.58 

2018/05/22 7.62 C 1.18 + 0.06 0.54 x 0.54 2.9 
J2 0.63 +0.03 0.83 0.54 x 0.54 

2018/11/18 8.65 C 1.09 + 0.05 0.61 x 0.19 6.7 
J2 0.53 +0.03 0.88 0.56 x 0.56 

2020/07/07 8.64 C 0.91 + 0.05 0.35 x 0.35 5.5 
J2 0.73 +0.04 0.88 0.56 x 0.56 

2021/07/07 8.64 C 0.92 + 0.05 0.37 x 0.37 4.7 
J2 0.63 +0.03 0.95 0.68 x 0.68 

2021/11/24 8.64 C 1.11 + 0.06 0.34 x 0.34 6.9 
J2 0.67 +0.03 0.95 0.69 x 0.69 

2022/07/05 8.64 C 1.04 + 0.05 0.27 x 0.27 10.2 
J2 0.59+0.03 0.97 0.67 x 0.67 

2022/09/06 8.64 C 1.09 + 0.05 0.39 x 0.13 15.8 
J2 0.52 + 0.03 1.00 0.60 x 0.60 


from 1.8 to 3.8 based on the values given in Marscher (1983); 
Vm is the turnover frequency in gigahertz (GHz); 0 is the angular 
size of the emitting region in mas; Sm is the flux density at the 
turnover frequency in Jy. This equation assumes a specific ge- 
ometry (e.g., a spherical emitting region) and that the magnetic 
field and particle energies are in equipartition. vm and S m are ob- 
tained from SSA spectrum model fit (Appendix D). @ and 6 are 
from VLBI model fitting (Appendix B). 

The uncertainty of the calculated Bssa can be expressed as 
OBssa = fBssa © Bssa. Assessing fgs, can be made through error 
propagation: 


fi = 5. Za) 4 (4.28) y 
Bssa Ym 6 Sm 


The magnetic field strength B calculated from Equation E.1 
is the value at the distance of the opaque core from the black 
hole. To provide some practical values at typical distances, one 
needs to convert a magnetic field strength measured at the loca- 
tion of the core (e.g., Bssa s$) to a value at 1 parsec from the black 


hole (Bıpc) usually involves understanding the scaling relation- 
ship of the magnetic field strength with distance in the jet. 

The magnetic field in AGN jets often follows a power-law 
decay with distance from the central engine, expressed as: 


h -Æ 
m-ai) 


where: B(h) is the magnetic field strength at distance h from the 
black hole, Bo is the magnetic field strength at a reference dis- 
tance ho (often the core or a specific emission region), £ is the 
power-law index that describes how the magnetic field decreases 
with distance. 

To convert Bssq (measured at the core or the region asso- 
ciated with the spectrum turnover) to Bj,-, one would need to 
know the distance of the core from the black hole in parsecs (or 
use the frequency to estimate this distance based on models of 
jet geometry) and apply the appropriate power-law index (£) that 
describes the magnetic field’s radial profile in the jet. 
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Table B.2. Parameters of the SSA fitting. 


Epoch Sin Vm Qthick Qthin Bssa Bipe Djer/M 
(Jy) (GHz) (mG) (G) (102G cm?/Mo) 
(1) (2) (3) (4) (5) (6) (7) (8) 
2009/10/15 3.70.1 8.3 += 0.3 1.4+0.3 -06+0.1 7.7+41.2 10.1 + 1.6 8.1413 
2009/10/16 3.6+0.1 8.2 +0.5 1440.3 -0.5+0.1 7.7423 10.1 + 3.0 8.1+42.4 
2009/11/15 4640.2 94+0.5 1.0+0.1 -12+0.2 6.4+ 1.9 8.3 + 2.4 6.6 + 1.9 
2010/08/19 3.3401 9.7405 0.80.2 -0.7+0.2 11.0+3.1 12.6 +3.5 10.14 2.8 
2011/01/15 3.30.1 8.2+03 0.80.2 —0.8 =+ 0.2 94+1.7 10.0 + 1.8 8.0 + 1.4 
2011/10/13 3.00.1 8&.1+0.6 0.80.2 —0.8 =+ 0.2 10.3 + 3.6 9.9 + 3.4 7.9 2.7 
2011/12/25 2.9401 7.9402 0.80.2 —0.8 + 0.1 10.9 + 1.6 10.1 1.5 8.11.2 
2017/04/10 2.5+0.1 59406 0.80.5 —0.6 + 0.2 7.2+3.8 3.8 + 2.0 3.0 1.6 
2017/10/01 2.3+0.1 65406 0.50.2 —0.8 + 0.3 10.8 + 5.4 5.5 = 2.7 4.4 + 2.2 
2017/12/18 2.3+0.1 62+0.8 1.9+0.5 —0.3 0.1 5.1 + 3.2 2.6 + 1.6 2.0 + 1.3 
2018/04/29 2.1401 5.0+0.7 19+0.6 -0.2+0.1 4.6 + 3.3 2.2+ 1.6 1.81.3 
2018/07/13 2.0401 520.6 0.50.2 -0.7 + 0.3 10.3 + 6.2 4.9 + 3.0 3.9 + 2.4 
2019/10/09 2.0+0.1 44403 0.40.3 -0.7 + 0.2 8.2 + 2.6 3.5 1.1 2.8 + 0.9 
2021/12/07 2.1401 65410 0.40.3 -0.8+0.5 10.6 + 8.6 3.7 + 3.0 3.0 + 2.4 
2022/11/05 2.3+0.1 12.2+2.1 0.20.1 -1941.2 15.8+13.5 4.9 + 4.2 3.9 + 3.4 
2023/11/07 2.6+0.1 82+04 100.6 -—0.8 =+ 0.3 14.1 + 3.9 13.3 + 3.7 10.7 + 2.9 
2024/02/20 2.4+0.1 84+09 0.50.1 -0.7+0.2 14.9+7.9 12.4 +6.5 9.9 =+ 5.2 


Note The uncertainties showed for fitting parameters including Sm, Vm, Qthick, aNd Qinin are fitting errors. The 
uncertainties of magnetic field feature measurements are transferred from the input parameters. 
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Fig. C.1. Temporal changes in the positions of jet components J1 and J2 in OH 471. J1’s positions are measured from the 2.3 GHz VLBI 
data, and J2’s positions are measured from the 8.4 GHz VLBI data. Proper motions for J1 and J2 are represented through fitted trajectories: 
My = —0.0025 + 0.0085 mas yr! (stationary), and uy = 0.040 + 0.003 mas yr! (4.4 +0.4c). 


The specific value of 6 can depend on the jet’s physical con- 
ditions and the theoretical model used. For example, a com- 
monly used model is the Blandford & Königl model, which as- 
sumes a conical jet structure with a magnetic field decreasing 
as h~! or h~?, depending on whether the magnetic field is domi- 
nated by the toroidal or poloidal component, respectively. In this 
study, 8 = | is assumed. 
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The Doppler factor 6 and redshift factor (1 + z) in equation 
E.1 play important roles in relating the observed properties to 
the intrinsic properties of the emitting region. The observed flux 
is boosted by 6**” and observed frequencies by ô, where æ is the 
spectral index. So the Doppler factor needs to be included in the 
magnetic field strength equation to relate the observed quantities 
to the intrinsic ones. Given an observed brightness temperature 
Tops and assumed Tint, 6 can be calculated as 6 = Tops /Tin. Under 
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the assumption of equipartition between particle and magnetic 
energy density, the upper limit of Tin can approximately to be 
5 x 10!° K (Readhead 1994). The (1 + z) factor accounts for 
converting between observed and rest-frame quantities. 

Converting the angular size observed at a certain frequency 
using VLBI to that at the turnover frequency involves under- 
standing the relationship between the observed frequency, the 
turnover frequency, and the angular size of the emitting region. 
The angular size of the emitting region at different frequencies 
can be related to the frequency by considering the optical depth 
(T). The angular size of the source is expected to decrease with 
increasing frequency in the optically thick regime due to the 
decreasing area of the emitting region that is above the optical 
depth threshold for significant emission. 

A common approach to estimate the angular size at the 
turnover frequency (O(vm)) from observations at another fre- 
quency (6(v)) involves the following equation: 


Am) = OY) a 


where vm is the turnover frequency, and v is the frequency of the 
VLBI observation, k is a constant that depends on the geometry 
of the emitting region and the spectral index. For a spherical, 
homogeneous synchrotron emitting source, k = 1/(1 + a) in the 
optically thick regime. 


Appendix F: magnetically arrested disk (MAD) 
F.1. Foundational Principles of MAD Theory 


The Magnetically Arrested Disk (MAD) theory proposes that 
strong magnetic fields built up near the black hole can arrest in- 
ward accretion, dominate the dynamics of the inner disk, and 
launch energetic jets. As matter accretes, it drags along and am- 
plifies the magnetic fields, leading to an accumulation of mag- 
netic flux that exceeds the disk’s diffusive capabilities. As a re- 
sult, the magnetic pressure becomes comparable to the ram pres- 
sure of the gas, significantly influencing the accretion flow and 
forming a highly magnetized, arrested disk state. Strong MAD 
magnetic fields promote relativistic jet launching by extracting 
rotational energy from the black hole. MAD disks can also ap- 
proach or exceed the Eddington luminosity limit, with magnetic 
fields playing a crucial role in collimating and driving very pow- 
erful jets. The process of magnetic flux accumulation allows 
episodic transitions to the MAD regime, characterized by the in- 
termittent production of jet components. 


F.2. Theoretical Framework and Equations 


The calculation of the poloidal magnetic flux, ®jet, that threads 
a parsec-scale jet is crucial in understanding the dynamics of 
jets near accreting supermassive black holes. This calculation is 
foundational for assessing the magnetic fields’ role in jet forma- 
tion, acceleration, and collimation. 

The poloidal magnetic flux, ®jet, is a measure of the total 
magnetic field passing through a certain area perpendicular to 
the jet’s flow direction. It is given by the integral of the magnetic 
field component (Bpo1), which is aligned with the jets length, 
over a cross-sectional area (A) of the jet: 


Diet = 1 Bpoi dA 
A 


(F.1) 


In the case of AGNs, it is often assumed that the mag- 
netic field has a significant poloidal component that can be de- 
scribed by models of magnetohydrodynamic (MHD) flows. The 
strength of the magnetic field (B) 1) at different distances from 
the black hole is estimated using observations and theoretical 
models. This can involve complex modeling of the jet’s emission 
properties, including Doppler boosting effects, and the physics 
of synchrotron radiation. The cross-sectional area (A) through 
which the magnetic flux is calculated needs to be modelled or es- 
timated based on observations. This often involves assumptions 
about the jet’s geometry (e.g., cylindrical or conical shapes) and 
size at the parsec scale. 


VLBI observations can resolve the jet structure and some- 
times directly measure the width of the jet at various points 
along its length. If the jet’s width can be directly resolved, this 
measurement can be used to estimate the cross-sectional area. 
The angular width of the jet, combined with the distance to the 
source, allows for a conversion to physical units (e.g., parsecs or 
light-years). For unresolved jets, the width might be estimated 
from the emission characteristics, such as the synchrotron self- 
absorption turnover, assuming a model for the jet’s emission 
(e.g., a conical jet model) and the geometry of the emitting re- 
gion. 


The cross-sectional area A depends on the assumed geometry 
of the jet. For a cylindrical jet, A = mr’, where r is the radius 
of the cylinder. For a conical jet, the area at a given distance 
from the central engine can be calculated if the opening angle 
of the cone and the distance along the jet axis are known, A = 
mh? tan?(6/2), where h is the distance from the jet apex, 0 is the 
observed angular size of the VLBI core. 


With estimates of Bpoi and the jet’s cross-sectional area, the 
poloidal magnetic flux can be calculated by integrating the mag- 
netic field strength over this area. 


The MAD condition can be quantitatively described by com- 
paring the magnetic flux Dgy threading the event horizon to the 
accretion rate M. A key dimensionless parameter that character- 
izes the MAD state is: 


E2) 


where: 

- Mpy is the magnetic flux through the horizon, 

- M is the mass accretion rate, 

- lg = GM /c? is the gravitational radius of the black hole (with 
G being the gravitational constant, M the black hole mass, and c 
the speed of light). 


In the MAD regime, the value of ¢ becomes large, indicating 
a strong magnetic flux compared to the accretion rate. The exact 
threshold value for ¢ indicating a transition to the MAD state 
can vary based on specific model assumptions. A typical value 
of ¢ = 50 is often used (Narayan et al. 2003; Zamaninasab et al. 
2014). 


While we cannot directly measure ®gy, we can observation- 
ally infer the poloidal magnetic flux threading parsec-scale jets, 
jet, by using the flux freezing approximation, i.e., Ogu ~ Dye. 
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F.3. Jet magnetic flux and magnetic field strength 


The magnetic flux jet is expressed in terms of the transverse- 
average magnetic field strength, B, as (Zdziarski et al. 2015): 


232 aryshB( + o)! 


Pie = 
Jet la 


Mpu (F.3) 


10°Mo 


= 8 x 10% f(ax)(1 + gel 


Ka 

1G 
This equation involves several parameters: 

- Ọja is the jet magnetic flux. 

-H =rg1 + - a’)!/?) represents the black hole event horizon 

radius. 

- s is a parameter related to the angular frequency of the field 

lines compared to the black hole’s angular frequency. 

- his the distance along the jet in parsecs. 

- lis the black hole angular frequency 

- a is the dimensionless black hole spin parameter. The quan- 

tity f(a*) = ry/(arg) is a function of the black hole spin. For a 

rapidly rotating black hole, f(a*) x 1. 

- Bis the magnetic field strength. 

-o = (6;/ s)’,s < 1 is the jet magnetization parameter, 

representing the ratio of Poynting flux to kinetic energy flux. 

(T@;) ~ 0.1, which led to similar o value with these in the lit- 

eratures (e.g., Chamani et al. 2021; Zamaninasab et al. 2014). 
To convert Bssa to B at a specific location (e.g., 1 parsec 

from the black hole), one would use a model that describes how 

the magnetic field strength varies with distance within the jet. 

This could be a theoretical model based on jet physics or an em- 

pirical relationship derived from observations (see discussion in 

Appendix D. Once Bıpc is known, it can be used in the calcu- 

lation of the jet magnetic flux using equation F.3. The magnetic 

flux would typically be calculated based on the cross-sectional 

area of the jet at 1 parsec and the magnetic field strength across 

this area. 
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